RMS Radio Source Contributions to the Microwave Sky 



S. P. Boughn and R. B. Partridge 
Department of Astronomy, Haverford College, Haverford, PA 19041; sboughn@haverford.edu 

t> ! ABSTRACT 

o 
o 



Cross-correlations of the WMAP full sky K, Ka, Q, V, and W band maps with the 
1.4 GHz NVSS source count map and the HEAO I A2 2-10 keV full sky X-ray flux 



2 , 

, map are used to constrain rms fluctuations due to unresolved microwave sources in the 

WMAP frequency range. In the Q band (40.7 GHz), a lower limit, taking account of 
fT) ; only those fluctuations correlated with the 1.4 GHz radio source counts and X-ray flux, 

corresponds to an rms Rayleigh-Jeans temperature of ~ 2^K for a solid angle of one 
i-£h ■ square degree. The correlated fluctuations at the other bands are consistent with a 

i i (3 = —2.1 ± 0.4 frequency spectrum. Using the rms fluctuations of the X-ray flux and 

2 ■ radio source counts, and the cross-correlation of these two quantities as a guide, the 

above lower limit leads to a plausible estimate of ~ 5/uK for Q-band rms fluctuations 
in one square degree. This value is similar to that implied by the excess, small angular 
scale fluctuations observed in the Q band by WMAP, and is consistent with estimates 

£> 1 made by extrapolating low-frquency source counts. 
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Subject headings: cosmology:observations — cosmic microwave background — radio continuim:galaxies- 
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^P. ' 1. Introduction 

Th e only available all-sk y survey for radio sources at frequencies above 5 GHz is the WMAP 



survey (jHinshaw et al. 20071 ). The WMAP team has cataloged several hundred radio sources de- 
tected at flux levels > 1 Jy in one or more of its five frequency bands; a refined search by Lopez- 
Caniego et al. (2007) has increased the number to nearly 400. Virtually all of these sources are 
already known from lower frequency radio catalogs such as FIRST, GB6, and the PMN survey in 
the southern hemisphere (see, for example, Trushkin 2003, Hinshaw et al. 2007, and Lopez-Caniego 
et al. 2007). The WMAP sources are a mix of flat spectrum sources, such as quasars, and the 
brightest "classical" radio galaxies that may have steep spectra, but retain adequate flux in one or 
more of the WMAP frequency bands. 

Since the WMAP sources do appear in lower frequency catalogs with high probability, it follows 
that correlating lower frequency catalogs with the WMAP images can either identify additional 
sources, as was found by Lopez-Caniego et al. (2007), or can reveal the statistical properties of 
sources too weak to be individually detected in the WMAP images. In particular, a cross-correlation 
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of the WMAP images with lower frequency catalogs can provide information on rms fluctuations 
induced by such sources at each of the WMAP frquencies, and thus average spectral indices for 
these subliminal sources can be deduced. Since the subliminal WMAP sources that dominate the 
fluctuations lie roughly in the decade 0.1 — 1 Jy, they are quite likely to appear in lower frequency 
catalogs, which generally reach to much fainter flux densities. This paper treats the correlation of 
WMAP images at 23 , 33, 41, 61, and 93 GHz with number counts from the NVSS 1.4 GHz survey 
(jCondon et al. 19981 ). Since many flat spectrum radio sources are known to be X-ray emitters, we 
also cross-correlate the WMAP images with the HEAO I A2 2-10 keV all-sky X-ray survey. 



2. Data Sets 

2.1. HEAO and NVSS Maps 

The versions of the HEAO and NVSS maps used in the present analysis were used previ- 
ously in a v ariety of analyses inc luding: an investigation of the large-scale structure of the X-ray 



background ( Boughn et al. 20 02); a determination of the large-scale bias of the X-ray background 



(Boughn fc Crittenden 2004ah : and measurements of the ISW effect (jBoughn fc Crittenden 20021 : 



Boughn fc Crittenden 2004b! ) . 



The NRAO VL A Sky Survey (NVSS ) is a flux limited survey (nominally at 2.5 mJy) at a 



frequency of 1.4 GHz ([Condon et al. 19981 ). It is complete for declinations above —40° and contains 
1.8 x 10 6 sources, with a mean source number density of 51.8 per square degree. The source counts 
were binned in 2 4,576 1.3° x 1.3° pixels. Th e pixelization we use is a quadrilateralized spherical 



cube projection (jWhite and Stemwedel 19921 ). Those pixels that were only partially contained in 
the survey region were omitted. There are systematic number count offsets that occur in several 
declination bands that coinc ide with discontinuous changes in the rms noise levels in the survey 



(jBoughn fc Crittenden 20021 ). To correct for this, random sources were added to or subtracted from 
each pixel to eliminate the band structure. The resulting map shows no declination dependent 
structure at a level of < 1%. For comparison, the Poisson noise per pixel is 11%. In order to 
exclude Galactic sources as well as nearby clusters of galaxies, the map was masked with the most 
agressive mask (KpO) used by the WMAP team. In addition, regions that contained pixels with 
source counts greater than 4<r (corresponding to > 43% of the mean number of sources per pixel) 
above the average count of 87.5 were masked. This procedure cleaned the map of 10 "objects" 
located more than 10° from the Galactic plane. Among these are the Orion Nebula and the nearby 
Virgo, Perseus, and Fornax clusters. That the moderately nearby, rich Coma cluster of galaxies 
was not one of the regions cut out indicates to us that this masking removed only Galactic and 
nearby extra-galactic sources. As is the case with the HEAO source masking, the effect of this 
extra cleaning was primarily to reduce the noise and did not significantly change the value of the 
correlation with the WMAP data. There is a small residual dipole in the NVSS map that is much 
too large to be that due to the Earth's motion with respect to the comoving frame of the radio 
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sources and, therefore, is undoubtedly due to some large scale systematic. This dipole moment was 
fit and removed from the map. 



The HEAOl A2 data set (jBoldt 19871 ) we employ here was constructed from the output of 
two medium energy detectors (MED) with different fields of view (3° x 3° and 3° x 1.5°) and 
two high energy detectors (HED3) with these same fields of view. These data were collected 
during the six month period beginning on day 322 of 1977. Counts from the four detectors were 
combined, then pixelized with the same projection used for the NVSS source counts. The combined 
map has a spectral bandpass (corre sponding to quantum efficiency > 50%) of approximately 3 — 



17 keV (jJahoda Mushotzkv 19891 ). but the counts have been converted to equivalent flux in the 
more standard 2 — 10 keV band. Because of the ecliptic longitude scan pattern of the HEAO 
satellite, sky coverage and therefore photon shot noise is not uniform. However, the variance of 
the cleaned, corrected map, 2.1 x 10" 2 (TOT cts s" 1 ) 2 , is significantly larger than the variance 
of pho ton shot noise, .8 x 10" 2 (TOT cts s" 1 ) 2 , where 1 TOT cts s" 1 ki 2.1 xlO" 11 erg s" 1 



cm -2 ( Allen et al. 1994 ). This implies that most of the variance in the X-ray map is due to 
"real" structure. For this reason, and to reduce contamination from any systematics that might be 
correlated with the scan pattern, we chose to weight the pixels equally. 

As with the NVSS map, the HEAO map was masked with the WMAP KpO mask and then 
more aggressively masked by removing all pixels within 20° of the Galactic plane and within 30° of 
the Galactic center. In addition, large regions (6.5° x 6.5°) centered on 92 nearby, discrete X-ray 



sources with 2 — 10 keV fluxes larger than 3 x 10 -11 erg s" 1 cm -2 ( Piccinotti et al. 19821 ) were 
removed from the maps. Around the sixteen brightest of these sources (with fluxes larger than 
1 x 10~ 10 erg s" 1 cm" 2 ), the masked regions were enlarged to 9° x 9°. We also used the ROSAT 



All-Sky Survey (RASS) Bright Source Catalog (jVoges et al. 1990 ) to identify additional bright 0.5- 



2 keV sources. While the RASS survey has somewhat less than full sky coverage (92%), it has a 
relatively low flux limit that corresponds to a 2 — 10 keV flux of ~ 2 x 10 -13 erg s" 1 cm" 2 for a 
photon spectral index of a = —2. Every source in the RASS catalog was assigned a 2 — 10 keV flux 
from its hard channel (0.5-2.0 keV) flux by assuming a spectral index of —3 < a < —1 as deduced 
from the hardness ratio within this band. There were 34 sources with fluxes exceeding 3 x 10 
erg s" 1 cm" 2 that were not in the Piccinotti source list, and these were also masked. The HEAO 
I A2 map itself was searched for "sources" that exceeded the nearby background by a specified 
amount. Thirty additional such sources at a level of ~ 3 x 10 -11 erg s" 1 cm" 2 were identified and 
masked. These final two cuts amount to only 6.5% of the sky and, in any case, do not significantly 
change the values of the correlations presented in Sj3j The primary reason for masking nearby, 
strong sources is to reduce the noise in the correlations discussed in The final masked map has 
48% sky converage. Finally, several systematic foreground features were fit and removed from the 
map: a linear time drift of detector sensitivity, high latitude Galactic emission, the dipole induced 
by the earth's motion with respect to the X-ray background, and emission from the plane of the 
local supercluster. These corrections are discussed in detail in Boughn (1999) and Boughn et al. 
(2002). 
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2.2. Microwave Maps 

The microwave data consists of the five frquency bands, K (centered at 22.7 GHz), Ka (33.0 
GHz), Q (40.7 GHz) V (60.6 GHz), and W (93.4 GHz), of the three year WMAP intensity maps 



(jHinshaw et al. 20 07). These maps were also converted from the available Res 9 HEALPix format 
to the same 1.3° x 1.3° quadrilateralized spherical cube pixelization of the other two maps. In 
order to exclude the strong microwave emission from the Galactic plane, all maps were masked with 
the agressive WMAP KpO mask. We note that, in addition, this mask removes all extragalactic 
microwave sources with fluxes greater than 1 Jy. For comparison, we also analyzed the maps without 
masking these sources, as discussed in §5.21 All five maps were corrected for the cosmic microwave 
background emission, which would otherwise add noise to the cross correlations, by subtracting 
the WMAP ILC CMB map. While this correction reduces the noise slightly, it does not have a 
significant effect on the lower limits deduced in §H The corrected maps still have a significant 
component of high latitude Galactic emission, especially on large angular scales. Since fluctuations 
due to extragalactic sources are predominently on small angular scales, the maps were high passed 
filtered by simultaneously fitting and removing the dipole, quadrupole, and octopole moments. A 
Galactic secant law and the Ha map compiled by Finkbeiner (2003), a template for Galactic free- 
free emssion, were included in the fits in order to remove additional Galactic emission. The resulting 
masked maps were then smoothed with a 24° FWHM Gaussian filter ignoring the contributions of 
masked pixels. These smoothed maps were subtracted from the corrected maps in order to provide 
additional high pass filtering. For unmasked pixels near large masked regions, the high pass filtering 
is somewhat reduced; however, only 1% of these pixels have the effective smoothing area reduced 
to a value below 0.4 times the nominal and no pixels had effective smoothing areas reduced below 
0.11. In any case, the effects of this high pass filter on small scale fluctuations were minimal, as 
discussed in £15.21 



3. Correlations 



Since the primary objective of the present analysis is to constrain the fluctuations of the 
equivalent brightness temperature due to microwave sources, it might seem reasonable to simply 
compute the rms fluctuations of the various WMAP maps. However, even with the masking, fil- 
tering, and corrections discussed in £j2j the rms fluctuations are completely dominated by Galactic 
emission. In order to isolate that component due to distant extragalactic sources, we cross-correlate 
the microwave maps with our two templates of extragalactic sources, the NVSS radio source count 
map and the HEAO X-ray flux map, which is dominated by distant AGN. The redshift distribu- 
tion of the signals represented by both of these maps are broad wi th a peak redshift of z ~ 1 
(jCondon 1984 : iDunlop Peacock 19901 : iBoughn Crittenden 2004al ) . A standard measure of the 
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correlation of two data sets is the cross-correlation function (CCF), which in this case is defined by 

CCF(e) = ^-^2(S l -S)(T j -f) (1) 

hi 

where the sum is over all pairs of pixels i,j separated by an angle 9, Si is the number of NVSS sources 
(or the X-ray intensity) in the ith pixel, S is the mean number of sources (X-ray intensity), Tj is 
the WMAP equivalent thermodynamic temperature of the ith pixel, T is the mean temperature, 
and Nq is the number of pairs of pixels separated by 9. In order to interpret the cross-correlations 
of the two templates with the microwave maps, it is necessary to determine the rms fluctuations in 
the template maps as well as the cross-correlation of these two maps. The auto-correlation function 
(ACF) is defined, using the same notation as in Equation [H as 

ACF{9) = ±- J2(Si - S)(S, - S). (2) 

Figure 1 and Figure 2 depict the ACFs of the NVSS source counts and X-ray flux. We have 
corrected the ACF(O) point of Figure 1 for Poisson noise by subtracting 87.5, the mean number 
of sources per pixel. With this co rrection, the signal at all an gles is an indication of the angular 



clustering of NVSS radio sources (jBoughn k, Crittenden 20021 ). The X-ray ACF at zero lag has 



been corrected for photon shot noise, but not for the Poisson noise in the distribution of sources 
since the resolution of the map, 3.04 degrees, is significantly larger than a pixel. However, a model 
of X-ray clustering indicates that the fraction of ACF(O) due to Poission noise is approximately 
77%, coincidentally the same as for the NVSS ACF(O). Figure 3 is the NVSS/HEAO CCF. The 
(highly correlated) errors in these ACFs and th e CCF have not been plotted; however, the signal t o 



noise values are large, i.e., S/N > 10 for 9 < 3° (jBoughn et al. 20021 : iBoughn fc Crittenden 2004al ) 



and have an insignificant effect on the analysis of this paper. 

Figure 4 is the CCF of the corrected NVSS map with the WMAP Q-band map. It is clear 
that only the small angular scale fluctuations of the two maps are correlated. Since the noise in the 
WMAP maps is not well characterized at the level of interest in our analysis, ~ 1/uK, we chose to 
estimate the noise by computing 240 CCFs of maps that were rotated with respect to each other. 
Five rotations were performed about each of 48 uniformly spaced axes separated by ~ 21 degrees. 
The resulting distribution of CCFs is, to a reasonable approximation, Gaussian and so the error 
bars that are plotted in Figure 4 have the usual 1 a interpretation. The errors so determined 
are accurate at the 10% level, which is sufficient for the arguments in §H As is clear from the 
figure, these errors are significantly correlated. This CCF is typical of the the CCFs at the other 
WMAP frequencies. Figure 5 shows the dependence of CCF(0) on frequency. The curve represents 
a power law microwave spectrum typical of free-free emission, T oc i/~ 2,1 , and it is clear that the 
data are consistent with such a model, although there is no clear detection at the highest (W-band) 
frequency. The curve is not a fit to the data but is forced to go through the highest signal to noise 
point at 33 GHz. Formal fits to these data are discussed in SjJJ Figure 6 and Figure 7 are the 
corresponding cross-correlation between Q-band and X-ray maps, and the CCF at zero lag for all 
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Fig. 1. — The auto-correlation function of the NVSS map corrected as discussed in the text. The 
units are (counts/pixel) 2 . 
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Fig. 2. — The auto-correlation function of the HEAO 1 A2 2-10 keV X-ray map corrected as 
discussed in the text. 1 TOT cts s" 1 = 3 x lO -11 erg s 1 cm 2 . 
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Fig. 3. — The cross-correlation function of the X-ray and NVSS maps of Figures 1 and 2 
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Fig. 4. — The cross-correlation function of WMAP Q band map with the NVSS map. Errors are 
determined from 240 rotated versions of the data and are highly correlated. 
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WMAP frquencies. It is clear that the detection of these correlations is 2a at best for the K, Ka, 
and Q bands, with no significant detection at the V and W bands. 

4. Constraints on Microwave Fluctuations 

In order to constrain fluctuations in the microwave sky we assume that the microwave emssion 
M in a given pixel consists of a signal aN proportional to the number TV of NVSS sources in the 
pixels plus a signal C that is statistically independent of N, i.e., M = aN + C. Then the zero lag 
value of the cross-correlation of M with N is given by 

CCF(0) = (MN) = a{N 2 ) (3) 

where () indicates an ensemble average, in our case to be approximated by the average over the 
maps. (N 2 ) is just the NVSS auto-correlation function, which we have measured. Solving for for the 
unknown constant of proportionality gives a = (MN)/(N 2 ). Finally the mean square microwave 
fluctuation is 

(MM) = a 2 (N 2 ) + (C 2 ) + 2a(NC). (4) 

Substituting for the constant a and noting that, by assumption, {NC) = 0, the rms microwave 
fluctuations are given by 

om = ^/WWJW)TW) > (MN)/y/(N2). (5) 

Therefore, a lower limit to the microwave fluctuations is given by the WMAP/NVSS CCF(O) 
divided by the square root of the NVSS ACF(O). The error in this lower limit is just the error in 
the CCF divided by the same factor, ^/(N 2 ). Table 1 lists these lower limits for all the WMAP 
frequencies along with those deduced in similar fashion from the WMAP/X-ray cross-correlation. 
No attempt is made to interpret the values of the CCF(#) for 9^0 since there is no detectable 
signal in the NVSS CCFs and the finite angular resolution (3.04 degrees FWHM) of the X-ray maps 
is responsible for most of the signal at 9 > 0. 

It should be noted that the ACF(O) used to deduce the values in Table 1 are those corrected 
for Poisson fluctutations for both the NVSS and X-ray data. It seems reasonable to do this since it 
is far more likely that the simple model of proportional emission applies to sources that cluster with 
each other than the assumption that the Poisson fluctuations of the two populations are correlated. 
Had we not corrected the ACF(0)s for Poisson fluctuations, then the lower limits in Table 1 would 
decrease by a factor of 2. The consequences of these alternatives will be discussed in Sj5j 

While the lower limits derived from the X-ray data are quite noisy, it is interesting that they 
are fortuitously similar to the values derived from the NVSS data. It is also of interest to estimate 
the lower limit of microwave fluctuations by considering the microwave emission that is associated 
with either NVSS or X-ray sources. It's straightforward to expand our simple model by expressing 
the microwave emission M as a combination of emssion proportional to N, the number of NVSS 
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Fig. 5.— CCF at zero lag of all 5 WMAP maps with the NVSS map. 
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Fig. 6. — The cross-correlation function of WMAP Q band map with the X-ray map. Errors are 
determined from 240 rotated versions of the data and are highly correlated. 



-10- 




Fig. 7. — CCF at zero lag of all 5 WMAP maps with the X-ray map. 

Table 1. Lower Limits on Micrwave Fluctuations Deduced from NVSS and HE AO 
Cross-correlations. The units are of rms Rayleigh-Jeans temperature. 



WMAP Band 


Secondary Map 


Lower Limits 


K Band (22.7 GHz) 


NVSS 


3.31 ± 0.78 




HEAO 


2.57 ± 1.61 




Combined 


3.24 ± 1.06 


Ka Band (32.95 GHz) 


NVSS 


1.77 ± 0.39 




HEAO 


1.29 ± 0.71 




Combined 


1.77 ± 0.41 


Q Band (40.70 GHz) 


NVSS 


0.92 ± 0.31 




HEAO 


0.88 ± 0.52 




Combined 


0.85 +0.45 -0.26 


V Band (60.60 GHz) 


NVSS 


0.45 ± 0.22 




HEAO 


0.34 ± 0.35 




Combined 


0.40 +0.44 -0.15 


W Band (93.44 GHz) 


NVSS 


0.04 ± 0.29 




HEAO 


0.27 ± 0.48 




Combined 


0.24 +0.63 -0.22 
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source counts, a component proportional to X-ray flux X, and a component uncorrelated with 
either N or X, i.e., M = aN + bX + C. Then the two CCFs are (MN) = a(N 2 ) + b(NX) and 
(MX) = a{NX) + b{X 2 ). Since we measure all quantities in these two expression except a and 
b, it is straightforward to solve for them and then evaluate the rms microwave fluctuations as in 
Equation [5J i.e., 

2 (MAQ 2 (X 2 ) + (MX)(N 2 ) - 2{MN)(MX)(NX) 2 
aM = (N*)(X*)-(NX) + {C) - (6) 

This combination does not constitute one that in any sense maximizes the signal to noise; rather, 
the purpose is to arrive at an estimate of the rms microwave fluctuations that are correlated with 
either the NVSS number counts or the HEAO X-ray emission. The lower limits to this estimate 
of the microwave fluctuations, ignoring the uncorrelated (C 2 ) term, are also listed in Table 1. 
Because the expression in Equation [6] is non-linear, the errors are non-Gaussian, especially for the 
lower signal to noise entries. They were determined from the 68% confidence limits deduced from 
the 240 CCFs determined from rotated maps. That the lower limits implied by the combined data 
are in agreement with the those determined separately from the NVSS and X-ray data is consistent 
with the relatively high degree of correlation of these two data sets. It should be noted that the 
NVSS and HEAO maps were not masked in the same way, so if there is any dependence of the 
ACFs and CCFs on masking then the interpretation of their combination in Equation 6 will be 
complicated. We repeated the analysis including only those pixels present in both the NVSS and 
HEAO maps. The values in Table 1 changed by only a small fraction of the error (~ 1/3<t) with 
the errors correspondingly larger because of the smaller number of pixels. 

Figure 8 shows the frequency dependence of the lower limit to microwave fluctuations deduced 
from the NVSS/WMAP CCF(O). The fluctuations have been expressed as equivalent Rayleigh- 
Jeans temperature in order to facilitate the comparsion with a power spectrum model. Because the 
corrected WMAP maps are still dominated by Galactic emission, these maps are highly correlated 
with each other. The result is that the errors depicted in Figure 8 are also significantly correlated. 
For example, the correlation coefficient of the K and Ka band errors is 0.79, of the K and Q band 
errors is 0.61, and of the Ka and Q band errors is 0.85. The solid curve is a minimum x 2 At to a 
free-free {(3 = —2.1) power spectrum and the dashed line is a fit to a typical synchrotron (/? = —2.8) 
spectrum. The reduced x 2 °f the former fit is 3.3/4 and 8.6/4 for the latter. A free-free spectrum 
is certainly consistent with the data while a synchrotron spectrum is marginally inconsistent with 
the data. (The reason that the synchrotron fit is less consistent than might appear from Figure 8 is 
because of the highly correlated errors.) One might expect that the model curves would be straight 
lines; however, due to the differences in the WMAP beam patterns, each data point corresponds 
to a slightly different smoothing that results in slightly different fluctuations than if all beams 
were the same angular scale. If the spectral index is included in the fit, then the minimum x 2 At 
yields (3 = —2.11 ± 0.41 with a reduced x 2 °f 3.25/3. The amplitude of the solid curve in Figure 
8 is 1.2(^/40.7 GHz) _2 ' 11 / uK. This value is well below the level at which systematic effects in the 
WMAP are characterized, which is the primary reason for determining the errors from the data 
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themselves rather than trying to model the noise of the maps. Another reason is that the NVSS 
map has significant systematics of its own, as discussed in 

It should be noted that the lower limits deduced above assumed that the correlations were 
due entirely to the joint clustering of NVSS and microwave sources. If instead we assume that the 
correlations are due to the fact that the NVSS and microwave sources are the same sources, then 
it would be appropriate to substitute the total, i.e., not Poisson corrected, ACF(O) in Equation [6j 
In this case the the lower limits are all reduced by a factor of 2.1. The consequences of this will 
be discussed in the next section. In order to facilitate comparisons with other observations and 
models we express the above results as the rms temperature fluctuations in a one square degree top 
hat beam. This value depends somewhat on the relative contributions of clustering and Poisson 
noise to the fluctuations. For example, if the Poisson process dominates, then the implied lower 
limit in one square degree is 1.9(z//40.7 GHz)~ 2,1 ^K while if the fluctuations are dominated by 1/9 
clustering, then the lower limit becomes 1.4(z//40.7 GHz)~ 2,1 /xK. If we take as a more reasonable 
assumption that ~ 25% of the fluctuations are due to the clustering of sources, comparable to the 
cases of both NVSS and X-ray fluctuations, then the limit becomes ~ 1.8(^/40.7 GHz) _21 ^K. 

5. Discussion 

The results presented in $3] are definitive measurements of the fluctuations in the microwave 
sky that are correlated either with counts of NVSS radio sources or with HEAO X-ray flux. Even 
so, they are only lower limits on the total fluctuations in the microwave foreground and the values 
were deduced under the assumption that correlations were due to joint clustering of microwave, 
NVSS, and X-ray sources. The power law fit to the lower limits indicates that the spectral index 
is more consistent with free-free or other flat spectrum than with synchrotron emission, which 
shouldn't be surprising since high frequency observations will naturally be biased toward flatter 
spectrum sources. 

5.1. Estimate of Total Fluctuations 

It is possible to make a (slightly) educated guess as to the total level of microwave fluctuations 
by appealing to the example of NVSS/X-ray correlated fluctuations. The level of X-ray fluctuations 
correlated with NVSS fluctuations can be deduced in the same manner as in Equation [U i.e., 
(X 2 ) > (X N) 2 1 (N 2 ) . The value of this quantity is 1.6 x 10~ 3 whereas the total X-ray fluctuation 
(see Figure 2) is (X 2 ) = 0.0123. Hence W (X 2 ) is a factor of 2.8 larger than the derived lower 
limit. The same factor applies to the case of a lower limit of NVSS fluctuations derived from 
X-ray fluctuations because Poisson fluctuations happen to account for the same fraction of total 
fluctuations in both the NVSS and X-ray cases. As was pointed out in §21 the NVSS and X-ray 
fluctuations have similar redshift distributions. IF microwave sources also have a similar redshift 
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Fig. 8. — Lower limits of the rms microwave flux in units of Rayleigh-Jeans temperature for all 
5 WMAP frequencies. The solid curve is the minimum x 2 fit to the data of a j3 = —2.1 free- 
free spectrum and the dashed line is the fit of a (5 = —2.8 synchrotron spectrum. (See $4] for an 
explanation of the slight curvature of both lines.) 
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distribution, IF microwave fluctuations have similar ratios of Poisson to clustering components, 
and IF the uncorrelated component ( C in Equation 0]) comprises the same fraction of the X-ray 
fluctuations as of the microwave fluctuations, then one would expect the same 2.8 correction factor 
to apply. This would imply that the total microwave fluctuations would be ~ 5 (u/40.7 GHz)~ 21 /iK. 
While plausible, there is no direct evidence that the three above conditions are true. Therefore, 
this conclusion should be considered only as a plausibilty argument. 

It was pointed out above that had total ACF(0)'s, including the Poisson contributions, been 
used to determine correlated microwave fluctutations, then the derived lower limits would be a 
factor of 2.1 less than those listed in Table 1. A similar analysis using the NVSS/X-ray correlation 
to determine a lower limit to the X-ray fluctuations would have have resulted in a lower value by 
this same factor. However, the factor by which the actual X-ray fluctuations exceed the correlated 
component would have increased by precisely the same factor. Therefore, while the lower limits 
depend on whether the correlated fluctuations are clustering or Poisson dominated, the rough 
estimate of total fluctuations does not. 



5.2. Potential Sources of Systematic Error 

All of the analysis described above was for WMAP maps that were masked for all regions that 
contained microwave sources that exceeded 1 Jy. We repeated the analysis not excluding these 
sources (but retaining the mask of the Galactic plane). The K and Ka band lower limits both 
increase by a factor of 1.4 while the Q band lower limit increases by a factor of 1.6. Given the 
limited signal to noise, these factor are consistent with each other, and our conclusions about the 
spectral index of fluctuations in the microwave sky remain unchanged. 

There may be some question as to whether the the high pass filtering discussed in £j2] above 
reduces the CCF(O) and hence the lower limits derived in this paper. As an estimate of the 
significance of this effect, the same procedure described in Sj2] was performed on the NVSS map, 
which was then cross-correlated with the unfiltered NVSS map. The result was a negligible (1.5%) 
decrease from the ACF(O) value of the original NVSS map. This should be an indication of the 
level of attenuation expected from the high pass filtered WMAP cross-correlated with the NVSS 
map. The high pass filtered NVSS map was also cross-correlated with the HEAO map in order 
to estimate the attenuation expected from the high pass filtered WMAP cross-correlated with the 
HEAO map. This attenuation was ~ 2%, again negligible. Since the resolution of the K and Ka 
bands (0.88 and 0.66 degrees) are not negligible compared to the pixel size (1.3 degrees), the former 
test might not be appropriate for the filtered WMAP/NVSS CCF. However, when a filtered X-ray 
map (which has a FWHM angular resolution of 3 degrees) was correlated with the NVSS map, the 
CCF(O) is attenuated by only 1.5% with respect to the unfiltered CCF(O), again a neglible amount. 
So we conclude that the high pass filtering has a negligible effect on the above results. 
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5.3. Comparison with Other Measurements and Predictions 

There is evidence of excess, small angular scale, fl uctuations in the three year WMAP Q 



band data (jHuffenberger et al. 20061 ; iHinshaw et al. 20071 ) . Under their assumption of uncorrelated 



(i.e., Poisson only) fluctuations, it is straightforward to translate this excess to equivalent Q-band 
Rayleigh-Jeans temperature fluctuations in a one square degree beam. The Hinshaw et al. (2007) 
estimate implies rms fluctuations of Trj ~ 4.5/iK while the Huffenberger et al. (2006) estimate 
corresponds to Trj ~ 4.0/iK. That these two estimates are consistent with our rough estimate 
above should be tempered by the assumptions required to arrive at the latter. To truly unravel 
the excess, small scale fluctuations, a sensitive, blind radio survey at the high frequencies normally 
used in CMB studies would be very helpful. In the meantime, our results suggest no unpleasant 
surprises in the WMAP frequency range. 

CMB measurements at smaller angular scales are more prone to contamination by foreground 
sources, since the power spectrum, £(£ + 1)CV, of unclustered sources is proportional to £ 2 . There 
is a substantial literature on the effect of radio sources on the microwave power spectra reported 
by the CBI, DASI and ACBAR groups, for instance. This is nicely summarized by Toffolatti et 
al. (2005). Toffolatti and his colleagues also attempt to model the expected contributions to high 
£ fluctuations from various classes of foreground sources. Unfortunately, these predictions can not 
be directly compared to our measurement because the models, like the microwave measurements, 
exclude radio sources down to a limit well below the ~ 1 Jy cut-off we employed. For instance, the 
CBI group (Readhead et al. 2004; Mason et al, 2003) removed all sources down to ~ 30 mJy (and 
many weaker sources as well) from their images before calculating the power spectrum. Therefore 
the measured (Readhead et al. 2004) and predicted (Toffolatti et al. 2005) levels of fluctuation lie 
roughly a factor ten lower than our value projected to £ = 1000 and a frequency of 30 GHz. 

On the other hand, we can compare our value of ~ 5 (j^/40.7 GHz)~ 2,1 /iK to values of rms 
fluctuation predicted by Toffolatti et al. (1998), since these predictions exclude only sources below 
1 Jy. Figure 6 of that paper shows that predicted level of fluctuations on a ~ 1 degree scale is 
~ 4.5/iK. When we scale from 44 GHz to 40.7 GHz using a (3 = —2.1 temperature spectrum, and 
correct from the Gaussian angular scale Toffolatti used to a 1 degree top hat beam, we obtain 
~ 5.3/uK, remarkably close to our estimate. Further, if we compare temperature fluctuations shown 
in his figures 5 and 6 for 30, 44, 70 and 100 GHz, we find that the implied spectral index on scales 
of ~ 1 degree is (3 ~ —1.9, close to our best-fit value, and inconsistent with a normal synchrotron 
spectrum. We conclude that we have found no evidence in the WMAP frequency range of a large 
contribution to the microwave sky not captured in Toffolatti's 1998 models. 

In the decade since 1998, refined models of source counts at high microwave frequencies have 
been published by, among others, Cleary et al. (2005), De Zotti et al. (2005), and Waldram et al. 
(2007). Unfortunately, these authors do not provide calculations of the expected fluctuation level 
in the microwave sky. We can, however, compare their source counts with those of Toffolatti et al. 
(1998) at the appropriate frequency and in the relevant range of flux density. In general, the more 
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recent source counts fall roughly a factor 1.5 — 2.0 below Toffolatti's. Thus the expected level of 
fluctuations would be 20 — 40% lower than Toffolatti's 1998 value, but still in agreement with our 
results. Obviously, the roughness of our estimate does not allow us to discriminate between these 
various models. What we can say is that, up to 90 GHz at least, there are no surprises in the form 
of unexpected populations of high frequency sources not present in lower frequency radio catalogs. 

We thank K. Jahoda for the original HEAO I A2 data, J. Pober who initiated an early analysis 
of residual microwave fluctuations, and J. Peebles for helpful discussions. We acknowledge the use of 
the Legacy Archive for Microwave Background Data Analysis (LAMBDA). Support for LAMBDA 
is provided by the NASA Office of Space Science. This work was supported in part by NSF grant 
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